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ABSTRACT 

We present a description of the public code XIM , a virtual X-ray observatory. XIM can be used to convert hy- 
drodynamic simulations of astrophysical objects, such as large scale structure, galaxy clusters, groups, galaxies, 
supernova remnants, and similar extended objects, into virtual X-ray observations for direct comparison with 
observations and for post-processing with standard X-ray analysis tools. By default, XIM simulates Chandra 
and the International X-ray Observatory IXO , but can accommodate any user-specified telescope parameters 
and instrument responses. Examples of XIM applications include virtual Chandra imaging of simulated X-ray 
cavities from AGN feedback in galaxy clusters, kinematic mapping of cluster velocity fields (e.g., due to merg- 
ers or AGN feedback), as well as detailed spectral modeling of multi-phase, multi-temperature spectra from 
space plasmas. 

Subject headings: radiation mechanisms: thermal, methods: numerical, telescopes, galaxies: clusters: general, 
X-rays: general 



1. INTRODUCTION 

As the peaks of the matter distribution in the cosmic web of 
structure formation, clusters hold a special place in cosmol- 
ogy. The ever-increasing power of computers has made direct 
hydrodynamic simulation of complex astrophysical systems a 
reality. We can now simulate large cosmological boxes that 
include massive clusters of galaxies and incorporate complex 
physical phenomena (such as star formation, the energy in- 
put from black holes, metallicity injection from supernovae, 
to name a few). 

X-ray observations of clusters have proven to be a rich and 
powerful tool to study large scale structure. The launch of 
the Chandra X-ray Observatory in particular has revealed a 
wealth of new information on cluster dynamics and at the 
same time, posed a number of new puzzles, such as the ques- 
tion of how clusters maintain their core temperature against 
catastrophic cooling. 

While hydro simulations can address many of these ques- 
tions, it is vital that the output data from such simulations can 
be compared with X-ray observations in the most direct and 
faithful way possible. Generally, galaxy cluster are simple to 
simulate compared to, for example, galaxies, where complex 
ISM physics and star formation seriously complicate model- 
ing. For that reason, it is much easier to produce realistic 
simulations of clusters that can be faithfully compared with 
observations. 

This comparison requires computer codes that can turn the 
output from a hydro simulation into a virtual observation with 
a given X-ray telescope — a virtual X-ray observatory. 

In this paper, we present such a virtual telescope, the IDL 
code XIM , which is now publically available at the web site 
|http : / /www . astro . wise . edu/^heinzs/XIM| 

XIM is best suited to process output from Eulerian (grid- 
based) hydro codes and its direct application is the virtual 
observation of thermal emission from cosmic gas. While it 
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was written to directly simulate Chandra and the next genera- 
tion X-ray observatory IXO , it can be adapted to simulate any 
arbitrary X-ray telescope for a set of user-supplied response 
matrix files. 

Since X-ray telescopes are imaging spectrographs (each 
photon is time- and energy-tagged, thus allowing a time- and 
energy-resolved image to be constructed), the output from a 
virtual X-ray observation is a spectral imaging cube (such as 
would be obtained from an integral field spectrograph in op- 
tical spectroscopy, for example). Figure Q] demonstrates the 
operation of XIM in calculating such a data volume. 




FIG. 1 . — Sketch of XIM process: Turn a 3D data cube of density, temper- 
ature, velocity, metallicity, and filling factor into a (2+l)D spectral data cube 
(shown is an IXO simulation of the iron Ko line region of a radio galaxy like 
Cygnus A, see fj4j . The left cube is 450 kpc on a side, while the right cube 
has the same spatial dimensions on the x- and z-axes and spans an energy 
range from 6.2 to 6.8 keV on the y-axis (into the paper). 



As the astronomical community considers the possibilities 
presented by IXO , an ultra-high resolution spectrograph with 
unprecedented collecting area, it is vital that we construct 
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detailed models of what we can expect in terms of spectral 
line signatures from complex systems like galaxy clusters and 
supernova remnants. Hydrodynamic simulations and subse- 
quent virtual observation with IXO present one of the best 
ways to explore the incredible wealth of information such data 
will offer. 

Given that no observatory with anywhere near as good a 
spectral resolution, let alone combined with the planned large 
area nd and high angular resolution, has ever successfully 
flown, experience from past X-ray missions is clearly limiting 
our imagination, and more sophisticated tools will be neces- 
sary to properly plan for the advent of IXO . 

At the same time, long Chandra observations of diffuse 
objects can benefit greatly from prior numerical simulation 
and it is generally recommended that observers perform ray- 
tracing simulations when proposing long observations. XIM 
can perform such simulations and presents a powerful tool in 
constructing mock Chandra observations of possible targets. 

This paper is organized as follows: In fj2] we will describe 
the computational details of how XIM converts input data into 
virtual X-ray observations, while $3]presents some examples 
of how XIM could be used and a brief discussion of some of 
the limitations inherent in XIM . Section H] presents a brief 
summary. 

2. IMPLEMENTATION AND ALGORITHMS 

In this section, we will describe the X-ray imaging pipeline 
XIM. XIM consistst of a suite of publically available IDL pro- 
grams that automate the creation of simulated X-ray data for 
a range of satellites, with a particular focus on simulations of 
Chandra and IXO observations. 

It takes input from hydro simulations of astrophysical ob- 
jects with a focus on galaxy cluster- and cosmological simula- 
tions. Combined with a catalog of hydro or MHD simulations, 
it can thus be easily made into a virtual observatory that can 
be scripted and operated remotely (e.g., via a web interface). 

In the following we describe in sequence the steps taken by 
XIM towards a virtual X-ray observation. 

2.1. Data input 

XIM is directly suited for grid-based hydro simulations. All 
input must be gridded on rectilinear grids. Adaptive mesh 
refinement (AMR) is only possible in the sense of staggered 
meshes, not yet in the sense of completely independent cell 
sizes. 

The main input variables for XIM are the particle density, 
the gas temperature, and the gas velocity in the form of three 
dimensional arrays. Additionally, XIM can accept arrays of 
metal abundance (relative to solar) and volume filling factor 
of the emitting gas (these variables can also be be passed as 
scalars, in which case they will be applied to the entire simu- 
lation). 

In addition, physical coordinate vectors (and optional tele- 
scope pointing, roll angle, and off-set parameters) are used to 
calculate the emitting volume and angular size of each vol- 
ume element of the input data, given an object red-shift and 
cosmology. By default, XIM uses concordance cosmologi- 
cal parameters of Ha = 70km/s, ft — 1.0, fl\ = 0.7, 

^matter = 0.3. 

The simulation output is calculated on a user specified en- 
ergy grid, which can be set to correspond one-to-one to the 
energy channel mapping of the instrument. This allows the 
user to easily create full resolution spectra as well as multi- 
color images, or low-resolution spectra for reduced computa- 
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FIG. 2. — Twodimensional image of a model grid at default XIM resolution 
for Chandra simulations, as function of gas temperature and photon energy, 
for solar abundance, as calculated from APEC . 

tional expense (high resolution spectra at IXO resolution can 
be very costly to run). A uniform turbulent velocity broad- 
ening can be applied to the spectrum to simulate sub-grid ef- 
fects. A uniform foreground hydrogen column density can be 
specified for photo-electric absorption. 

2.2. Spectral models 
2.2.1. Thermal spectral calculation 

By default, XIM uses the publically available APEC model 
to calculate spectra from thermal plasma emission (assumin g 
gas in coronal equilibrium) for each cell (Smith et al. 2001a). 
APEC self-consistently calculates the equilibrium ionization 
balance for a thermal pla sma. Atomic data are taken from 
the ATOMDB using APED (ISmith et al.ll2001bl) and combined 
with bremsstrahlung continuum for all species. APEC output 
is provided as a table of spectral models. 

Given the fact that APEC itself provides a table of model 
spectra, and given that the computational expense of extract- 
ing a spectrum from APEC is high, the most economical 
method of calculating a large number of spectra with XIM 
is by creating an oversampled table of model spectra for the 
entire temperature range spanned by the simulation. Thus, the 
spectral projection in XIM is based on interpolation on a table 
of model spectra , on a logarithmic temperature and photon 
energy grid (see 92.31 . 

The default temperature binning of the table is 66 bins per 
decade in temperature (significantly oversampled with respect 
to the APEC model output and suffi cient for high-accuracy 
spectral interpolation, again see Ej2 . 31 > - 

Two sets of spectral tables are created: One for Helium and 
Hydrogen, assuming primordial cosmic Helium abundance, 
and one for heavier elements (with relative abundance of met- 
als fixed to the solar ratios). The metal abundance Z can be 
specified for each cell. A model grid at default Chandra res- 
olution is shown in Fig.|2] 

2.2.2. Other spectral models 

XIM is not limited to the implementation of APEC de- 
scribed above for spectral modeling. Any user-defined spec- 
tral model can be applied (e.g., simple powerlaw spectra) and 
an arbitrary set of spectral parameters and input data can be 
passed. 

However, because XIM is optimized for thermal emission, 
the spectral model must be such that the emission is propor- 
tional to some power of a density variable (this could be rel- 
ativistic electron pressure, which enters synchrotron emission 
as the 7/4 power) and a temperature variable (which controls 
the spectral shape and normalization, e.g., synchrotron age). 
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The output spectra must be provided in the form of a two- 
dimensional table in temperature and photon energy with log- 
arithmic spacing on the energy axis 

2.2.3. Foreground absorption 

XIM allows the specification of the neutral Galactic hydro- 
gen column density to calculate photo-electric foreground ab- 
sorption following the Wiscon sin Absorption Model (WABS 
iMorrison & McCammonll 19831) . 

2.3. Spectral projection 

The spectral projection along the line of sight is performed 
along an arbitrary viewing angle. For each cell of the input 
data, a unique spectrum is calculated for the cell temperature, 
cell density, and its proper Doppler shift, given the line-of- 
sight velocity of that cell. 

To project along an arbitrary line-of-sight, the input data 
arrays are rotated (using tri-linear interpolation) onto a new 
grid that orients the x-axis along the line-of-sight vector. 
The input velocity arrays are also projected onto the line- 
of-sight vector to calculate a line-of-sight velocity for every 
cell. Line-of-sight integration then simply proceeds along the 
x-axis of the new array. 

2.3. 1. Doppler shift and energy gridding 

The Doppler frequency correction from the emitting to the 
observed frame is 



Vohs — f^cmit — t'cmit -1 / " — ; T ( 1 ) 

y 1 + v-los/c 

If the frequency grid is logarithmically spaced, Doppler shift 
is a simple addition. This suggests that spectral calculations 
should be performed on logarithmic grids, which is the model 
adopted by XIM . On a logarithmic energy grid, Doppler shift 
is a simple linear shift along the grid, which is computation- 
ally trivial. 

Thus, given an output energy grid (either user-defined or 
fixed to the instrument channel map), XIM creates an over- 
sampled logarithmically spaced energy grid that is padded on 
both sides by the maximum Doppler shift encountered in the 
input data and takes into account the energy resolution of the 
instrument, additionally extending the grid to enclose a mini- 
mum of 99% of the energy at the grid edge. 

The energy resolution of the logarithmic grid is set to over- 
sample the target output grid by at least a factor of 2. Since 
the energy grid used by XIM is logarithmic, most of the en- 
ergy range is oversampled by a larger factor. 

Doppler shift is performed by linear interpolation between 
the nearest two integer shifts along the energy axis that strad- 
dle the actual Doppler shift value of each cell. This accurately 
reproduces the mean line centroid energy and leads to a small 
error in line width that is well below the energy bin width and 
well below the intstrument resolution. 

The accuracy of this procedure is demonstrated in Fig. [3] 
which shows the line centroid energy determined by XSPEC 
for a gaussian spectral line with varying Doppler shift virtu- 
ally observed with IXO (to sufficiently high signal-to-noise 
to eliminate statistical error). 

2.3.2. Line-of-sight integration 

In addition to the Doppler shift correction, the line-of-sight 
integration must interpolate the spectral model for each cell 
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FIG. 3. — Accuracy of the numerical Doppler shift procedure relative to 
the actual Doppler shift. Shown are XSPEC fits to the line centroid (dots) 
for a range of Doppler shifts applied to the gaussian line and then virtually 
observed with IXO (actual Doppler shifted centroid energy shown as solid 
line). The centroid is determined to very high accuracy. Also shown are the 
IXO energy resolution of 2.5 eV and energy channel width of 0.5 eV. 

from the table provided by the spectral model calculation de- 
scribed in $2.2\ 

This is done by linear interpolation between the two near- 
est temperature spectra in the model table. Given the de- 
fault dense sampling of the spectral model grid (66 temper- 
ature bins per temperature decade), this method provides a 
smoothly varying, accurate description of the spectra as a 
function of cell temperature. 

Given both temperature and Doppler calculations, the line- 
of-sight integration is then performed on a cell-by-cell basis. 
For an input grid of size N x x N y x N z and a logarithmic 
energy grid with N\ e bins, this involves N x N y N z N\ c calcu- 
lations, which can be very time consuming and lends itself 
to trivial parallelization. A simple implementation of paral- 
lelization of this step is included in XIM . 

The line-of-sight integration is performed before re- 
gridding to the instrument pixel scale to preserve the maxi- 
mum amount of velocity detail from the original hydro simu- 
lation. However, for high-resolution spectra, as encountered 
with IXO , memory requirements can become extreme. In 
some cases, it is therefore advisable to re-sample the input 
grid to lower resolution, appropriate for the lower angular res- 
olution of the telescope (however, this procedure will lead to 
an under-estimate of the Doppler width of the lines). 

To keep memory requirements at a minimum, spectra are 
accumulated on a pixel-by-pixel basis (i.e., the line-of-sight 
integration is performed along the line of sight first, then 
across the image). Thus, only the input data and the spectral 
cube of dimension N x N y N\ e have to reside in memory. 

2.4. Virtual observation 

The output from the line-of-sight projection encapsulates 
the object's cosmological red-shift, but it is still instrument- 
indepent (though often the energy grid will be chosen to be 
suitable to a particular instrument to keep computational ex- 
pense and memory requirements at a reasonable level; e.g., 
for virtual Chandra observations, emphasis should be placed 
on spatial resolution, while for virtual IXO obsevations, the 
energy grid should be at the highest resolution possible). 

Starting from the output file of the line-of-sight inegration, 
virtual observations with different target telescopes and in- 
struments are performed in a number of steps described be- 
low. 

2.4. 1 . PSF simulation 
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The line-of-sight integrated spectral cube is convolved 
with the telescope point spread function (PSF). XIM allows 
user specified PSF models. Note that off-axis effects on the 
PSF and vignetting are currently not included in XIM , i.e., it 
is assumed that the PSF is uniform across the field-of-view. 

The default PSF used for IXO simulations was extracted 
from the IXO simulator simx, which provides a monochro- 
matic PSF (specified at 6keV). 

The default PSF for Chandra simulations was extracted 
from MARX ray-tracing simulations and is interpolated in en- 
ergy with a native spacing of 0.5 keV. The low-energy resolu- 
tion is 0.5 arcseconds (see Chandra users guide for more in- 
formation on the energy dependence of the PSF). More accu- 
rate Chandra imaging simulations, including vignetting and 
off-axis effects, can be ach ieved using the optional MARX in- 
terface described in ^2.4.51 

2.4.2. Convolution with instrument responses 

A critical component of realistic X-ray simulations is the 
incorporation of the effective area of the telescope as a func- 
tion of photon energy (encapsulating mirror effective area and 
detector quantum efficiency, typically encoded in the ancillary 
response file in X-ray astronomy) and the redistribution func- 
tion of photon energy (reflecting the spectral resolution of the 
telescope and detector, typcially encoded in the response ma- 
trix file). This happens in three steps: 

• Given a response file (following the CAL/GEN/92-002 
FITS specifications), XIM either re-grids the telescope 
response to the logarithmic energy bin on the data 
energy axis and the user-supplied output energy grid 
on the the detector axis, or, if XIM uses native tele- 
scope energy channels, the line-of-sight projection is 
re-gridded from logarithmic onto detector channels. 

• Next, the line-of-sight projected, PSF-convolved data 
cube is re-gridded onto physical detector pixels for the 
specified telescope/detector. 

• Finally, XIM convolves the spectral cube with the in- 
strument response onto the specified output energy grid, 
resulting in a spectral cube of count rates per energy bin 
per pixel as detected by the instrument. 

2.4.3. Virtual detection 

XIM assumes Poisson noise for counting statistics and 
produces an IDL data file that includes the detector count 
rate (without statistical error) and a counts file for the user- 
specified exposure time that includes Poisson errors. 

In addition to IDL output, XIM can produce a FITS events 
file that can be analyzed with standard X-ray astronomy tools 
like f tools and ds9 . 

XIM can also produce a single FITS spectrum for the en- 
tire observation. Post-processing with the spectral genera- 
tor of XIM allows the application of IDL mask arrays to se- 
lect sub-regions of the spectrum (or select different weighting 
for different pixels for the combined spectrum) and produce 
FITS spectra that can be analyzed with XSPEC , isis , and 
sherpa . It is possible to create multiple different spectral 
files for different regions from the same ximulation (e.g., an- 
nuli around a cluster center for deprojection analysis). 

2.4.4. Backgrounds 



XIM adds sky background and instrumental background to 
the virtual observation. These can be user-specified spectra or 
standard background estimates included for Chandra and the 
IXO calorimeter XMS. 

Sky background is treated like source photons and is folded 
through the instrument response (though assumed to be uni- 
form across the field-of-view, so no PSF convolution is per- 
formed). Instrument background is added directly to the out- 
put counts without convolution with the response. 

The default background for Chandra was extracted from 
the black sky background file for post-2000 data. The IXO 
background estimate is based on the most recent estimates 
available on the telescope wiki page and will be continually 
updated as more sophisticated estimates become available. 

2.4.5. Virtual Chandra observations with MARX 

XIM comes with default Chandra observing parameters 
(the latest aimpoint response files should be downloaded by 
the user from the Chandra web archive). This alows high- 
fidelity simulations of objects near the aimpoint. 

However, for observations of objects far off-axis, PSF ef- 
fects and vignetting should be taken into account. Thus, for 
the most realistic virtual Chandra observations, XIM uses 
MARX (the official Chandra simulation software) to perform 
ray-tracing simulations of the observation. The output is again 
a fully compliant Chandra FITS events file that can be an- 
alyzed using ciao and the Chandra pipeline like regular 
Chandra data. 

3. DISCUSSION 

Having briefly described the steps taken by XIM to convert 
a grid of data taken from a hydro simulation into a virtual X- 
ray observation, we will present some examples of possible 
applications (the details of which will be presented in a sepa- 
rater paper) and discuss some of the limitations. 

3.1. Examples of applications 
3.1.1. Chandra imaging of radio galaxies in clusters 

Before going into more detail, we list several examples of 
possible applications of XIM for Chandra : 

• A good demonstration of the power of X-ray imag- 
ing is the unsharp-m asked image of the Perseus cluster 
dFabian et al.1 12003) that shows a series of concentric 
ripples (identified as sound waves and a weak shock). 
While direct one-to-one simulations of a Perseus are 
not possible given the difficulty in reproducing the de- 
tails in cluster weather that likely lead to the asymmet- 
ric appearance of the cluster, it is nonetheless very use- 
full to examine hydro simulations for whether they can 
reproduce the observed surface brightness fluctuations 
induced by these waves. 

• X-ray astronomy makes us e of object colors fo r selec- 
tion effects. For example, (iForman et al.ll2007h used a 
color cut on a long Chandra observation of the Virgo 
cluster to detect a weak shock driven into the cluster by 
the radio galaxy M87 at the center. 

• As a further example, a cosmological simulation con- 
taining multiple galaxy clusters can be examined for 
the optical set of bands within which to select galaxy 
clusters, depending on the instrument used. 
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Multi-color imaging is particularly easy and fast using XIM 
, given that only a small number of spectral channels have 
to be calculated (APEC calculates the integrated photon lumi- 
nosity for a given energy band, so no interpolation is involved 
when calculating broad band images). 

Fig.|4]shows the simulated Chandra surface brightness and 
the unsharp-masked image for a hydro simulation of Cygnus 
A. The details of the h ydrodynamical simulation are described 
in iHeinz et all (120061) and we limit this discussion to a very 
brief summary of the features of th is simulation: We used the 
publically available FLASH code dFrvxell et al.ll2000l) to sim- 
ulate the injection of a super-sonic, powerf ul jet into a galaxy 
cluster (taken from (Springel et alJ l2Q0lT) '). The jet power 
(10 46 ergs s _1 ) and the cluster temperature (~ 5keV) and 
central density match the canonical parameters for Cygnus A 
very well (e.g., citepwilson:02). When placing the simulated 
cluster at the red-shift z = 0.0561 of Cygnus A (picking a 
time step that reproduces the observed angular sizes of the X- 
ray cavities and the shock in Cygnus A) the simulated Chan- 
dra surface brightness is in good agreement with the actual 
Chandra observation (to within about 10%). 

The power and importance of virtual X-ray imaging can be 
demonstrated from the unsharp-masked image: The appear- 
ance of ripples in the wake of the outgoing shock wave is clear 
from this image. It implies that X-ray surface brightness rip- 
ples can be excited in the wake of one single outgoing shock 
wave as a result of jet dynamics and do not necessarily require 
multiple episodes of activity of the central radio galaxies. 

Given the abundant data on X-ray cavities in clusters in 
the Chandra archive, a statistical treatment of bubble prop- 
erties has now become possible. When computing the ob- 
servability and imaging properties of a sample of bubble ei- 
ther from simulations or analytic models, it is critical to em- 
ploy faithful renderin g of the virtual X-ray data. Using XIM , 
Brueg gen et al.l (120091) have demonstrated the importance of 
line-of-sight projection for estimating bubble position and 
size from observations. 

3.1.2. High-resolution IXO spectroscopy of a virtual clone of 
Cygnus A 

Using the same hydro simulation described in 33.1.11 we 
used XIM to simulate a high-resolution spectrum of a power- 
ful radio galaxy like Cygnus A observed with IXO . A more 
detailed description of the power of high-resolution spec- 
troscopy for the study of radio galaxies will be published in a 
companion paper. Here we wish to demonstrate the power of 
such observations using a simple example. 

The high spectral resolution of IXO will allow detailed 
kinematic mapping of the expanding cavities of radio galax- 
ies. Figure|5]shows three spectra, one through the cluster cen- 
ter and two through the expanding cavities. The line emission 
is due to the K-a lines from FeXXV and FeXXVI (Helium- 
and Hydrogen-like iron), which emit abundantly in hot clus- 
ter atmospheres like that of Cygnus A. 

The two spectra taken through the cavity centers show clear 
evidence for Doppler-shifted line emission from the front and 
back walls of the cavity and it is possible to simply read off 
the expansion velocity from the spectrum to be 614 km s -1 
from the separation of the three right-most peaks of the cyan 
spectrum. These three peaks are all due to the resonance line 
of FeXXV from the rear wall (left peak at 6.341 keV), the 
front wall (right peak at 6.367 keV) and the fore- and back- 
ground cluster emission (middle peak at 6.354 keV). The ex- 
pansion line-of-sight velocity in the simulation is u s hoii ~ 




FIG. 4. — Simulated Chandra image of Cygnus A (top) and unsharp- 
masked version of the same image (bottom). Three colors correspond to 0.5- 
2.0 keV (red), 2.0-5.0 keV (green), and 5.0-10.0 keV (blue). The image x and 
y dimensions are 440 x 330 kpc. The unsharp-masked image clearly shows 
ripples reminiscent of those seen in, for example, Perseus A. Note that the 
images shown are calculated with infinite signal-to-noise to demonstrate the 
appearance of ripples in X-ray images behind shock waves. Given a proper 
exposure time, it would be trivial to determine the detectability of features at 
a given surface brightness. 

650 km s -1 . 

More sophisticated model fitting would allow us to deter- 
mine the emission measure distribution at different velocities, 
possibly allowing a deconvolution of the internal dynamics 
of the shell (which is not moving at uniform velocity). This 
could be used to, for example, determine the rate at which gas 
is flowing back along the shell as well as detecting the cre- 
ating of high entropy gas around shells as direct evidence of 
heating and feedback. 

Other obvious examples of using high-resolution X-ray 
spectroscopy in the context of galaxy clusters include detailed 
APEC model fits to determine the temperature and metallicity 
distribution of the cluster gas. Demonstrations of the power of 
this method are beyond the scope of this paper, but the XMM- 
Newton obervations of galaxy clusters that gave rise to the 
second cooling flow problem (ultimately leading to the cur- 
rent paradigm of AGN heatin g of clusters) provi de evidence 
of the power of this method dPeterson et alj |2003) even at the 
lower resolution and throughput of XMM-Newton . 

It is also possible to incorporate and test multi-phase mod- 
els versus observations. This can be easily achieved by super- 
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FIG. 5. — Simulated IXO spectra of Cygnus A for three lines of sight 
(green/yellow: through the left cavity; black: Through the cluster center; 
cyan: Through the right cavity). The three most prominent peaks of the cyan 
spectrum show the approaching and receding cavity walls, as well as the clus- 
ter fore- and background emission. The offset relative to the cluster emission 
(black) is due to some rotation in the cluster and the jet inclination relative to 
the line of sight. The image inset shows the simulated image of the cluster. 

imposing multiple XIM simulations (one for each phase). 

3.2. Limitations 

Finally, we will briefly discuss some of the current limi- 
tations on the software that are not likely to be remedied or 
implemented in the near future. 

3.2.1. Grid geometries 

Based on the projection method, the code can only handle 
input data on regular grids that can be described by a sin- 
gle, rectungular three dimensional matrix. This excludes true 
AMR grids (though it allows for staggered grids) and, natu- 
rally, Lagrangian formulations from direct application. To use 
XIM with AMR or SPH output data, they will have to be re- 
gridded onto a regular grid. The same is true for non-cartesian 
grids. The examples shown above were all calculated with the 
AMR code FLASH and then re-gridded onto a regular carte- 
sian grid. 

3.2.2. Line-of-sight projections 

Direct integration along the coordinate axes is computation- 
ally simpler and faster, which is the method employed by XIM 
. This implies that all input arrays have to be rotated prior to 
the spectral projection. As with regular projections through a 
rectangular volume for arbitrary line-of-sight vectors, miss- 
ing values lead to the appearance of projection effects (e.g., 
the volume is rendered as a projected box. Currently, XIM 
does not allow for period boxes to be rendered without data 
loss. It is possible to clip the data, however, to avoid the ap- 
pearance of projection artifacts from missing values. 



3.2.3. Spectral models 

As already described, XIM incorporates the APEC thermal 
plasma emission code for plasma in coronal equilibrium. The 
method used to interpolate spectra assumes that the spectral 
shape depends only on one parameter (the temperature or an 
equivalent variable that can be passed instead of temperature) 
and that the emissivity is proportional to the square of the 
density (or to some power of an equivalent variable that can 
be passed instead) 

XIM allows user-defined spectral models as long as they can 
be written in such a way that they depend on two parameters 
(one that determines the spectral shape, such as temperature 
and one that affects the normalization through a powerlaw de- 
pendences). As long as this is the case, it is straight forward 
to incorporate other models. For example, synchrotron emis- 
sion could be included by having the spectral shape (index 
and/or cutoff) determined by one parameter (which is passed 
as the variable t, but could in fact be the spectral age of the 
plasma) and the overall normalization is determined by the 
non-thermal pressure (which would be passed as the parame- 
ter n). Spectral shifts would be incorporated through the line- 
of-sight velocity v. 

Clearly, this method is somewhat limiting in the way that 
spectral models can be formulated. For example, photo- 
ionized plasmas will be difficult to implement. 

3.2.4. Imaging restrictions: Vignetting and off-axis 

As already mentioned above, XIM assumes a uniform PSF 
across the field-of-view and neglects vignetting. For IXO , 
these extent of these effects is currently not well known. For 
Chandra , the best way to incorporate these effects is to use 
the MARX option to pipe the spectral cube through the tele- 
scope ray tracing software. This will produce the most realis- 
tic Chandra simulations. 

4. CONCLUSIONS 

We have presented a description of the public code XIM . 
The code produces virtual X-ray observation from hydro sim- 
ulation of astrophysical plasmas. The code is best suited for 
the simulation of thermal emission, using the APEC model. 
XIM incorporates detailed telescope parameters and responses 
for Chandra and IXO . Examples of the application for 
XIM include multi-color imaging, detailed surface-brightness 
mapping (e.g., using unsharp-masking), high-resolution kine- 
matic studies of the dynamics of galaxy clusters, groups, and 
supernova remnants, and multi-temperature, multi-phase high 
resolution spectroscopy of space plasmas. Current limitations 
include the restriction to regular, rectangular grids, spectral 
models that can be parameterized in a fashion equivalent to 
thermal emission, and restriction to on-axis imaging perfor- 
mance. 
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